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Abstract 

Domains of condensed-phase monolayers of chiral molecules exhibit a variety 
of interesting nonequilibrium structures when formed via pressurization. To 
model these domain patterns, we add a complex field describing the tilt de- 
gree of freedom to an (anisotropic) complex-phase-field solidification model. 
The resulting formalism allows for the inclusion of (in general, non-reflection 
symmetric) interactions between the tilt, the solid-liquid interface, and the 
bond orientation. Simulations demonstrate the ability of the model to exhibit 
spiral dendritic growth. 
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I. INTRODUCTION 



Several experiments performed over the last decade have explored the myriad complex 
domain shapes formed during the growth of condensed-phase Langmuir monolayers [|TJ. 
An interesting case is the pressure-induced growth of condensed-phases whose constituent 
molecules are chiral. Here the growth shapes range from simple spiral dendrites to intricate 
chiral fractal structures [@-§|]. More recently, experiments with racemic monolayers have 
exhibited initially achiral dense-branched growth, followed by the addition of chiral "hooks" 
at the edges of the initial domains |J. 

In the same decade, the invention of the phase-field model provided a means for studying 
not only these isolated growth elements (such as dendrites), but also for elucidating the 
character of the global morphology of these domain growths. In a standard phase-field 
model of solidification, the phase of the material is represented by a real order parameter 
field whose time evolution is coupled to a diffusive field 0-0. As the liquid solidifies, the 
latent heat released by the advancing solid front acts as a source to the diffusion field. In the 
sharp- interface limit, this formulation recovers the traditional free-boundary solidification 
equations. Simulations of this model are able to reproduce the dense-branched structures 
seen in experiments with isotropic materials. In addition, anisotropic terms may be added 
to the standard phase-field model allowing for simulation of the stable dendrites seen in 
experiments with anisotropic materials [|T0l -[T2]|. Still, the chiral dendrites seen in experiments 



remain beyond the grasp of these models. 



In a recent advance |L3| , the standard phase-field model was generalized to include infor- 
mation regarding the orientational order of the "solid" phase (for example, the orientation 
of the crystal axes of a solid) by incorporating this information directly into the phase of 
a complex order parameter. Prior to this work, the implicit "crystal axes" were fixed with 
respect to the coordinate axes. While the introduction of a complex order parameter al- 
lowed the "crystal axes" to bend dynamically, chiral structures were still excluded by this 



model. Even with the inclusion of non-reflection symmetric anisotropy (a prerequisite for 
chiral growth) the dendrites failed to spiral. This anisotropy merely lead to the growth of 
non-reflection symmetric dendrites growing at constant velocity |T4| . 

Some clues as to the physical origin of spiral dendrites have been provided by equilib- 
rium studies of textures in Langmuir monolayers (single layers of amphiphillic molecules at 



an air- water interface) fllSyiBj . In these studies, many of the condensed phases are observed 
to have hexatic orientational order. In one, called the "LS-phase" , the average direction 
of the hydrophobic tails is perpendicular to the air-water interface. In others, the average 
tail alignment has a component in the plane of the air-water interface. Landau free energy 
descriptions of the interactions between the tilt and the bond orientation have been formu- 
lated which successfully describe the phase transitions among the various condensed phases 
IB 



With this motivation, our strategy will be to add a complex order parameter representing 
the tilt to an anisotropic complex-phase-field model, and to couple the tilt direction both to 
the bond orientation and to the orientation of the solid-liquid interface. As we shall see, the 
resulting model is indeed capable of supporting the growth of spiral dendritic structures. We 
interpret this as proof that these patterns can be understood as being the result of the incon- 
sistency of the tilt preferring to point both along a particular direction with respect to the 
interface normal and along a particular direction with respect to the crystalline anisotropy. 
That is, we conclude that spiral dendrites can arise from frustration of the tilt field. 

II. THE MODEL 

A complex-phase-field model describing diffusion-limited solidification is defined by the 
equations [[L3|]: 

d* _ 6F[*,U] 

dt 5^* 1 1 
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F = F^, = J dr k 2 |W| 2 + |^| 2 -2|^| 4 + |^| 6 + T|^| 2 tanh(Aff) (3) 

where ^ = |^|e lAre is an N-fold rotationally symmetric complex order parameter whose 
argument indicates the orientation of the crystal axes with respect to the x-axis. The first 
term in equation (||) is an energy cost for spatial inhomogeneity, which gives rise to a surface 
tension proportional to k. The remainder of equation (H) represents a "double-well" whose 
minima at =0 (liquid) and |\&| = 1 (solid) are tilted by coupling to the diffusion field U 
so as to favor the appropriate phase. 

Physically, the bond orientational order leads to an anisotropy in the surface tension. Pre- 



vious workers have included these anisotropic effects in various ways p. 10, 13). Kobayashi 



18| and Wheeler, Murray and Schaefer [|19j include anisotropy by allowing the surface ten- 



sion coefficient to depend on the local orientation of the gradient of the order parameter. 
McFadden, et. al. later provided a thermodynamic formulation of this method, as well 



as an asymptotic analysis of the sharp- interface limit. We adopt this method of including 
anisotropy, defining k as a function of both the crystal axis orientation and the normal to 
the interface. For an N-fold anisotropy we define 

«(*) = k {l + 77 K |Uf |V|#|T cosN(6 n -6)}, (4) 

with the normal to the interface (pointing from the solid into the liquid) defined, 9 n = 
tan -1 (Syl^/I 2 /^!^/! 2 ). The resulting dendrites will grow along the N-fold crystal axes. 

To describe the "tilt" of the molecules, we define a new complex order parameter 
$ = |$|e"^ whose magnitude reflects the local degree of tail alignment and whose argu- 
ment indicates the direction of tilt in the plane of the air- water interface. We again assume 
the time-dependence 

rf$ _ 5F[V, U] 



Experimentally, only liquid-condensed phases exhibit tilt. To incorporate this fact, consider 
the energy 



1 mag f]mag 



|rff(|$| 2 - |^| 2 ) 2 (6) 



which for large values of r] mag effectively locks |$| to \^\ (interpreting |$| = as the non-tilted 
phase and |$| = 1 as the tilted phase). In addition, the tilted phases are known to exhibit 
direction-locking at various angles with respect to the bond orientation. This experimental 
information may be incorporated via an energy 

F dir = r] dir J df \<5>\ 2 \V\ 2 {1 -cos N{<j)- 6 -9 dir )}, (7) 

where Q^ T is the equilibrium angle between the tilt direction and the crystal axis orientation. 

Note that the formulation just described is not meant to capture the full physics of 
the phase transitions among the various tilted phases. Previous workers have formulated 
Landau theories describing the phase transitions among the various condensed phases (those 
possessing orientational order) ||i6| , |T7| . While one could attempt to generalize this theory to 
describe transitions among the various condensed phases (both tilted and non-tilted) and the 
expanded phases (those with no orientational order and hence no "crystal axes"), we are not 
interested in transitions amongst all these phases. The above formulation is only meant to 
capture the physics of the transition from a liquid-expanded ("LE") phase with = to a 
single specific tilted condensed phase in which = |<3>| = 1. In the context of the literature 
surrounding tilted hexatic phases for achiral molecules, Odir — corresponds to an "L 2 "-hke 
tilted phase (where the tilt points to the nearest neighboring molecule), Odir = to an 
"L^'-like tilted phase (where the tilt points to the next nearest neighbor). Any other value 



of dir is akin to an "L^'-like phase |H| which has dynamic reflection-symmetry breaking. 
If the molecules are chiral, there will be no chiral symmetry breaking transition and Odir will 
take some specific value. 
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Finally, we add an interaction between the tilt angle at the interface and the interfacial 
orientation; this is done via the inclusion of the term 



In the sharp-interface limit, this interaction becomes a boundary condition on the angle of 
the tilt , 0, at the interface. To see this, consider a small, nearly planar portion of the 
interface. Since U varies most rapidly across the interface, we may approximate U as a 
function of only the coordinate normal to the interface, v = 2 cos a + ysina (where a is 
the angle from the x — axis to the gradient of U). Then d z U = \ (d x — id y ) U ~ ^e~ ta d v ll, 
and Fjj<$, becomes proportional to cos (</> — (3 — a). So we expect the tilt to make an angle 
(5 with respect to the interface normal at equilibrium. For the growth process, we expect 
this condition to be most strongly enforced at the dendrite tip as the diffusion gradient is 
enhanced there. Note that a nonzero value of (3 can only arise for achiral molecules. 

We have used the phase-field model just described as a computational stand-in for the 
"true" sharp- interface model (as defined by experimental observations). In the limit of a 
sharp interface, we expect the magnitudes of the order parameters \^\ and |$| to be locked 
to their equilibrium values. That is, there is a region of liquid-expanded phase in which 
l^l = |$| = surrounding a region of tilted liquid-condensed ("solid") phase in which 

= |$| = 1, with a physically "sharp" interface separating the two regions. Meanwhile, 
the phase angles of the order parameters (the bond direction 9 and the tilt direction 0) are 
free to vary continuously in the bulk "solid". The relationship between the tilt direction 
and the bond direction has been observed experimentally, and is approximated here by the 
energies @ and (|7J). While the boundary condition on the bond orientation cannot (as 
yet) be observed experimentally, boundary conditions on the tilt at the interface have been 
observed |T5|] . We approximate these tilt boundary conditions with the energy (|). 
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III. SIMULATIONS 



We have performed numerical simulations of the above model (defined by equations 
([]]), (0), and (||) ) with a free energy 

, $,U) = J dr 4| V$| 2 + + F mag + F dir + F U9 (9) 

where the surface tension due to $ is the N-fold symmetric function given by equation 
(H) and the surface tension due to <£> is isotropic. For convenience we rescaled the model 
to have length units of 4 = DT/k, (the "kinetics-limited diffusion length"), time units of 
11/ D (the related diffusion time), and an interface thickness of e (see Ref. |2(J for details of 
rescalings). Our present computing resources limited our simulations to 200x230 gridpoints 
(on a triangular lattice), which does not sufficiently resolve the interface. We used a three- 
fold bond orientation order parameter (N = 3) to help minimize competition between the 
implicit six-fold lattice anisotropy and the explicit anisotropy of the bond orientation (for 



further discussion, see Ref. |T3"|). 

Clearly, we do not pretend to present quantitatively accurate simulations of these spi- 
ral growths. Our goal here is to check that the set of interactions described above can 
indeed give rise to the spiral dendrites seen in the experiments. Further progress would 
require more fully resolved simulations and systematic procedures for reducing the effects of 
lattice anisotropy. We refer the concerned reader to the important work of Boesch, Mueller- 



Krumbhaar, and Shochet |21j which presents a detailed prescription for controlling the effects 
of lattice anisotropy in computer simulations of phase-field models. We view such careful 
studies as a necessary postlude to the phenomenological model presented here and as a 
prerequisite for any quantitative comparisons with experiment. 

In general, the simulations below used e = 0.4, A = 10, T = 0.6, N = 3, T$ = 0.75r, 
k$ = 10~ 4 K ,r] mag = 1.0, and rjdir = 1.0. Initial seeds were circles with radii of 6 with 
uniform ^-distributions whose three-fold directions are obtained by rotations of ^ from the 



+?/-axis; the initial (^-distributions are described below. The undercooling at infinity was 
fixed at a value A (i.e. U (oo) = —A) which varied in the simulations below. To prevent 
overcrowding in the figures, tilt vectors are shown for approximately one tenth of the points 
for which |\&| > 0.50 (the "solid"), resulting in domains which are sometimes slightly more 
planar than the fully sampled interfaces. Bond orientations made an angle 0<&r with the tilt 
(to within a tenth of a degree) and do not appear in the figures. 

The three-armed spiral dendrite of Figure la was achieved with a strong anisotropy 
(compared to the relatively slow dendrite growth speed, A = 0.50 and rj K = 2.0), a chiral 
boundary condition on the tilt direction at the interface (riu® = 2.5, (3 = ^j) and a tilt 
locked parallel to the local crystal axis orientation (6^ = 0, analogous to the "L2" phase in 
the literature fL5f). Within the first 100 timesteps, the initially radial 0-distribution relaxed 
to the domain structure shown in the inset of Figure 1, which is the three-fold analogy of the 



tilt domains seen experimentally in droplet textures of hexatic condensed phases ||15| . In this 
case, the chirality of the material was expressed in the boundary condition on the tilt which 
meant that the tilt rested at an angle of ~ clockwise from the normal to the interface. This 
resulted in a clockwise bending of the crystal axes via the energy i^ ir and hence a clockwise 
turning of the growth direction. 

As expected, there is a parameter regime in which the anisotropy is too weak (for a given 
dendrite growth speed) to support stable dendrites, resulting in a dense-branched structure. 
The combination of a higher dendritic growth speed (A = 0.90) and a lower anisotropy 
strength (j] K = 0.25) produced the pattern shown in Figure lb (all other parameters are 
identical to those in Figure la). This pattern resembles at least to some extent the chiral 
dense-branched structures seen in bacterial colony growth f22 |. 

An interesting scenario is the case where the tilted phase is chiral while the constituent 
molecules themselves are achiral. [For some discussion of chiral symmetry breaking by achiral 
molecules in Langmuir monolayers and Langmuir-Blodgett films see References [|T^|,^3J .] This 
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case is peculiar in that a given "L^'-like phase can produce either a right-handed or a 
left-handed spiral, depending on the orientation of the domain walls of the early tilt field 
with respect to the crystal axes. An example of this is shown in Figure 2 where stable 
dendrites (A = 0.65, i] K = 2.0) grow under the influence of an achiral tilt boundary condition 
(t)u$ — 2.5, P — 0) in an "Ll"'-like tilted phase (9 dir = =f). Constant angular displacements 
were added to otherwise radial ^-distributions to obtain the tilt domains (formed during the 
first 100 timesteps) shown in the insets. In Figure 2a, the tilt in the upper domain is locked 
j radians counterclockwise from the crystal axis along the +y-axis. As the dendrite tip 
emerges along the +y-axis, the tilt rotates clockwise to follow it. To remain in the "L'^'-like 
phase, the crystal axes (and hence, the growth direction) of the added material also turn 
clockwise. In Figure 2b, the tilt in the upper domain is locked j counterclockwise from 
the crystal axis oriented below the +x-axis. This time, as the dendrite tip emerges from 
the upper domain the tilt must rotate counterclockwise to follow it. The growth direction 
will then turn counterclockwise as the crystal axes rotate to remain in the existing "L^'-like 
phase. 



IV. CONCLUSION 

In summary, we have introduced a complex order parameter representing the tilt degree 
of freedom in an anisotropic complex-phase-field solidification model, and added terms to the 
free energy which represent interactions between the tilt, the solid-liquid interface, and the 
bond orientation. Though these added terms have not been derived from thermodynamic 
identities, we appeal to the sharp interface model suggested by the experimental observations 
of tilted hexatic domains as a motivation for these terms. 

We have shown that our model allows for the simulation of various chiral scenarios which 
produce spiral dendrites and chiral dense-branched structures. We note in passing that we 
have had to allow for significant flexibility in the tilt-field and bond orientation in the bulk 
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of the material. While this is reasonable for liquid crystal growths, presumably most other 
crystalline systems (e.g. metallic) would be much more rigid and would probably preclude the 
possibility of spiral growth. Nonetheless, this demonstration of a set of sufficient conditions 
for the formation of these structures should help further the understanding of the intricate 
global morphologies seen in experiments on diffusion-limited growth over the last decade. 
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FIGURE CAPTIONS 



FIG. 1. Growth shapes for chiral tilt boundary conditions (r]u^ = 2.5, (3 = = f) in an "L2"-like 
tilted phase {0^ ir = 0) . (a) Stable three-fold spiral dendrites for A = 0.50, r\ K = 2.0. 
(b) Chiral dense-branching for A = 0.90 and rj K = 0.25. 

FIG. 2. Stable three-fold spiral dendrites for normal boundary condition on tilt (r)u<f> = 2.5,(3 = 
0) in an "LI"' -like tilted phase (9 dir = =f), with A = 0.65, r] K = 2.0. Early tilt domain 
formations are shown in the insets. 
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